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Abstract. The dynamics of a Bose-Einstein condensate is studied theoretically in a combined periodic plus 
harmonic external potential. Different dynamical regimes of stable and unstable collective dipole and Bloch 
oscillations are analysed in terms of a quantum mechanical pendulum model. Nonlinear interactions are 
shown to counteract quantum- mechanical dephasing and lead to phase- coherent, superfluid transport. 

PACS. 03.75.Lm Tunnelling, Josephson effect, Bose-Einstein condensates in periodic potentials, solitons, 
vortices and topological excitations - 05. 45. -a Nonlinear dynamics and chaos 



1 Introduction 

The study of transport properties of ultra-cold atoms in 
corrugated potentials has become an intensely discussed 
topic since the first experiments with Bose-Einstein con- 
densates (BECs) in optical lattices almost a decade ago 
[T1|^1|51|11|51|S1|7||5 1|^ The observed or pre- 

dicted phenomena are often discussed by concepts bor- 
rowed from the phenomenology of extended systems. This 
approach led to the characterisation of superfluid and in- 
sulating phases and phase transitions |lHT3j , modulational 
instability P"^IITSj . and dissipative behaviour How- 
ever, different dynamical regimes like small- amplitude os- 
cillations, dephasing instabilities, and Bloch oscillations 
require different models and analogies for their explana- 
tion. 

In this work we approach the problem from a differ- 
ent, somewhat holistic point of view and treat the BEC 
in the external potential as a finite dynamical system. By 
mapping this problem onto a simple pendulum model, we 
are able to explain different dynamical regimes as well as 
stabilisation and destabilisation mechanisms in a unified 
approach. Specifically, we consider a cloud of ultra-cold 
bosonic atoms in a one-dimensional (ID) optical lattice 
in a classical tight-binding approximation with additional 
harmonic trapping in the lattice axis. We further confine 
the analysis to the situation of a sufficiently large num- 
ber of atoms per site that quantum fluctuations may be 
neglected . Typical experiments with hundreds of atoms 
on the central site certainly satisfy this condition but 
probably the theory still remains valid with much lower 
atom numbers. This system has been discussed before 
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in many experimental and theoretical works. A dynam- 
ical "superfluid-insulator" transition has been predicted 
in Ref. 5 for a BEC on the basis of a modulational insta- 
bility caused by nonlinear interactions. A similar effect of 
insulating behaviour, however, was observed in noninter- 
acting fermions The latter was interpreted as a very 
different mechanism in terms of a semiclassical pendulum 
model. Experiments with BECs j^f l^ and further numer- 
ical analysis [TTHT^ gave ambiguous results in showing 
reduced mobility of bosons without revealing the mecha- 
nism. 

In this paper we study the BEC in a combined har- 
monic and lattice trap by exploiting a mapping of this 
system to a quantum mechanical pendulum model. This 
exact mapping of the lattice dynamics of the noninteract- 
ing system to a simple quantum pendulum model estab- 
lishes two effects: A separatrix in the semiclassical phase 
space leads to two separate regions of qualitatively differ- 
ent dynamics (see Fig. Q] and discussion in Sec. [2|). Fur- 
thermore, dephasing occurs due to quantum-mechanical 
wave packet motion. Adding non-linear interactions intro- 
duces two additional effects: Far away from the separatrix, 
the nonlinearity counteracts the quantum dephasing of the 
linear problem and thus stabilises wave packet motion. We 
understand the emerging coherent wave-packet dynamics 
as a signature of superfluid transport in the sense dis- 
cussed in the recent literature (see, e.g. Refs. |Hll?)lll()llllj V 
Close to the separatrix and depending on the strength of 
the nonlinearity, a dynamical instability destroys coher- 
ent wave packet motion. We thus obtain a unified view 
of such different phenomena as dynamical instability, de- 
phasing, and superfluidity. While a dynamical instability 
has been predicted by a different mechanism in Ref. 5l 
and the effect of the separatrix has been discussed in con- 
juction with fermions in Ref. we believe that quantum 
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dephasing and its suppression by nonlinear interactions 
has not been discussed before 19 . In addition to the sta- 
bilisation of dipole oscillations, we also predict a so far 
undescribed regime of coherent wave packet motion above 
the separatrix. This regime could be exploited to gener- 
ate the recently proposed atomic gas at negative kinetic 
temperatures |20] . 

2 Non-interacting atoms 

We begin with dynamics of non-interacting atoms, gov- 
erned by the Schrodinger equation with the following 
single-particle Hamiltonian 

In Eq. (QJ M is the atomic mass, Vq the depth of the op- 
tical potential, d the lattice period, and uj the frequency 
of the harmonic confinement. If the lattice potential is 
large compared to the recoil energy Er = 2h^T:'^ / d?M ^ 
we may use the tight-binding ansatz for the Schrodinger 
equation, ilj{x^t) = ai{t)ipi{x)^ where ipi{x) are the lo- 
calised Wannier states. This leads to a system of coupled 
linear equations for the complex amplitude a^(t), 

ihai = ^I'^ai - ^ (a/+i + ai-i) , (2) 

where u = Mco^cP and J is the hopping matrix element, 
uniquely defined by the depth of the optical lattice s = 
Vo/Er as J/Er - 53/4exp(-2y^)/V^. 

A particularly transparent description of the dynamics 
governed by Eq. Q is obtained by mapping it to the math- 
ematical pendulum Indeed, introducing the function 
(/)(6>,t) = l/V^T.iailt)exp{iW), the system of Eq. 
reduces to the Schrodinger equation for the quantum pen- 
dulum with the Hamiltonian 

H = ^L^-Jcos{0), L = -i—. (3) 

This problem is related to the Mathieu equation p!^24ll25l 
which is solved by well-known special functions [TTj . 
The full advantage of the representation however, is 
the easily accessible interpretation in terms of pendulum 
dynamics. A characteristic feature of the classical pen- 
dulum is the existence of a particular trajectory - the 
separatrix, which separates the vibrational and rotational 
regimes of the pendulum, see Fig. C] The notion of the 
separatrix can be well extended into the quantum problem 
[2H1- It is associated with the critical angular momentum, 
or the critical site index of the original problem, 

r = 2{J/vf^^ . (4) 

The existence of a critical /* has been indicated in labora- 
tory experiments 6 , where the authors excite the sys- 
tem by suddenly shifting the harmonic trap by a dis- 
tance Ax. Then, for /q = Ax/d < I* (for the lattice 
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Fig. 1. Phase space of the classical pendulum in the variables 
/ and 0. The dash-dotted lines show the separatrix for /* = 
40, separating rotations from oscillations. These correspond 
to Bloch and dipole oscillations of a BEC, respectively. Thin 
(grey) lines indicate classical trajectories with initial conditions 
^0 = and lo = 13, 27, 50, —50. The quantum time evolution 
of a Gaussian wave packet is schematised by the 50% contour 
line of the corresponding Wigner function j21j . After initial 
displacement from equilibrium position to lo = 27, = 0, 
solutions of Eqs. are shown in anti-clockwise order for the 
quantum pendulum at ^ = (thick green contours) and the 
interacting case at ^/J = 1 (thin red contours) at times as 
indicated. 



parameters used in the cited experiment T = 134) the 
wave packet oscillates around the trap origin, while for 
lo > r it sticks to one side of the parabolic poten- 
tial and the centre-of-mass position can never reach the 
equilibrium position. For ^ = these dynamical regimes 
are illustrated in Fig. [21 The characteristic frequency of 
the wave-packet oscillation is given by the pendulum fre- 
quency f2{l) 1221 . Below the separatrix (/ <C /*) we have 
f2{l) « i?o = {i^Jf^'^/h = cj(M/M*)i/2, where M* is 
the effective mass of an atom in the lowest Bloch band. 
At the separatrix i7(/*) = and above (/ /*) we have 
f2{l) ^ vl/h. Note that the dynamics of the atoms for 
/ > r can be viewed as Bloch oscillations of the atoms in 
a (local) static field F = vl^/d with the Bloch frequency 
Qbo = dF/h. 

In addition to the effect of the separatrix one can also 
see the effect of dephasing in Fig. O which smears out 
the oscillations of the wave packet as time goes on. This 
can be related to the non-equidistant spectrum of the 
quantum pendulum, which is inherited from the nonlin- 
ear frequency dependence f2{l) of the classical pendulum. 
If r 1, a short-time description of the dephasing can be 
obtained by solving the equations of motion of the clas- 
sical pendulum for an ensemble of trajectories with ini- 
tial conditions scattered over the phase volume ~ 27rfieff 
with heff = (i//J)^/^ = 2/r. As a result we obtain a t^- 
exponential decay for the oscillations of the mean coor- 
dinate and momentum of the atoms. The dephasing is 
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Fig. 2. Dynamics of non-interacting atoms = 0). Gray scale 
images show the time evolution of the squared wave func- 
tion after the ground state of the system has been suddenly 
displaced by the distance Ax at t = 0. The left and right 
columns show real space and momentum space plots, respec- 
tively, for Ax/d = 8 (dipole oscillations) in the upper row and 
for Ax/d — 24 (Bloch oscillations) in the lower row. The lattice 
parameters correspond to J = 2.4 • 10~^Er, v — 3.2 • 10~^Er 
(hence, /* = 17). The time axis is scaled by the period 27r/i7o 
of small-amplitude pendulum oscillations. Momentum is scaled 
by the reciprocal lattice constant kh — 2Ti/d. The multiple 
peak structure of the momentum-space plots is due to the pres- 
ence of a periodic potential. 



conveniently quantified in terms of the quantity 



which has been introduced as an "order parameter" in 
Ref. 5 . In fact, due to the normalisation condition 
"^iCLidi = 1, we find ^ ^ 1 when the site-to-site phase 
fluctuations are small and ^ ^ in the presence of strong 
phase fluctuations. The upper left panel of Fig. 01 shows 
the decay of ^ during dipole oscillations due to dephasing. 

Concluding this section we note that the calculations 
were done by using the continuous nonlinear Schrodinger 
equation with a lattice depth of s = 12.16. For this depth 
of the optical potential and the considered initial displace- 
ment, the results obtained within the tight-binding ap- 
proach (not shown) practically coincide with the depicted 
ones. The deviation between the solutions appears only 
for the initial shift Iq larger than Imax ^ A/u, where A 
is the energy gap between the ground and first excited 
Bloch bands. If Iq > Imax the Landau-Zener tunnelling 
takes place and the single-particle dynamics of the atoms 
is a superposition of the Bloch and dipole oscillations. 



Fig. 3. The same as in Fig. |5|but for finite nonlinearity g 
1.55 • 10"^^;^. 



3 Interacting atoms 

We shall analyse the case of interacting atoms in the frame 
of the ID Gross-Pitaevskii equation, 

.f^di^ = HiP{x,t) + g,Dmx,t)\^i^{x,t) , (6) 

where giD ^ ashoo±N^ as is the s-wave scattering length, 
u;± the radial frequency and N the total number of atoms. 
The tight-binding version of © reads as 

ihai = ^I'^ai - ^ (a^+i + a/_i) g\aifai , (7) 

where g = gm J \ tlJi{x)\^dx ~ ashuj±N/d. 

The main result we want to report in this work is that a 
weak nonlinearity can compensate the dephasing and the 
wave packet follows the classical trajectory of the pendu- 
lum without dispersion. This is illustrated in Fig. [HI and 
the right column of Fig. ^ 



3.1 Variational approach 

In order to estimate the amount of nonlinearity required 
to convert the quantum dynamics of the pendulum into 
the 'classical dynamics', we use the Gaussian variational 
ansatz of Ref. ^Q^. For the quantum pendulum this 
amounts to a semiclassical approximation that can ac- 
count for dephasing [SI]. In this approach, the wave packet 
is parametrised as 



ai 



{t) VAexp 



(8) 
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Fig. 4. Decay of the order parameter ^ over time for different 
initial displacements /o below /* = 40 (upper panels) and above 
r (lower panels). The panels on the left hand side correspond 
to the noninteracting case {g — 0) and the panels on the right 
to g/J 1. 



where A is a normalisation constant. Then the centre of 
the wave packet L{t)^ the dispersion 7(t), the velocity ^(t), 
and the dephasing parameter 5{t) satisfy Hamilton's equa- 
tions for the effective Hamiltonian 

i/e. = |(i^ + ^)-Jcos^e- + ^, (9) 

where rj = 1/27^ + 7^(5^/8 and the pairs of canonical vari- 
ables are {L^6) and (7^/8,(5), respectively. Thus we have 

hL = — J sm6e~^ ^ W = vL ^ hj = JjS cosOe~^ ^ 

h6 =JcosO (^-sA e-^ + -^^^ - V . (10) 
V7 / V^7 

The non-dispersive dynamics of the wave packet depicted 
in the lower row of Fig. [21 implies the (quasi)periodic dy- 
namics of the variables L, 7 and J. In fact, for the cer- 
tain range of the nonlinearity g and harmonic confinement 
there is a stable periodic orbit in the four-dimensional 
phase space of the system which comes through the 
point (5 = 0. The condition for the existence of this peri- 
odic orbit is approximately given by the condition 

glv = ^fl2, (11) 

which means that the last two terms in the equation for (5 
cancel each other. Examples of the discussed stable peri- 
odic orbits are shown in in the upper two panels of Fig. [SI 
In lower panels of Fig. Owe plot the stability regions of the 
orbit together with the estimate (|TT|) . The bright regions 
correspond to the quasiperiodic dynamics with small devi- 
ations of b and 7. In the grey (red) regions the deviations 
are large and in the dark (blue) regions the dephasing 
J increases without bounds. It is worth noting that the 
variational ansatz (|S|) becomes invalid as soon as the or- 
bit is unbounded or badly bounded. On the other hand, 



if (5(t) is captured around (5 = and 7 is not too small as 
it occurs near the centre of the stability island, we have 
exp[— 7^(t)] 1 and Eq. (|Tn|) reduces to the equation of 
motion for a classical pendulum. 

Let us estimate the minimum strength of nonlinearity 
needed in order to suppress dephasing. A coarse estimate 
may be obtained by the requirement that the equilibrium 
width of the wave packet is compatible with Eq. ([TT|) . 
The minimum value ^min can be found by requiring that 
the width 7 obtained from Eq. (|TT|) is equal to the non- 
interacting equilibrium width from setting (5 = and ^ = 
in Eq. (|10|) . We obtain the condition 




dephasing 6 dephasing 6 




Fig. 5. Upper panels: Poincare cross section of the effective 
system Q for /* = 40 in the oscillating regime with /o = /*/2 
(left) and the rotating regime with /o = 4/* (right). The peri- 
odic orbit is located in the centre of the stability island. Lower 
panels: The stability region of the depicted periodic orbit in 
the (7,^)-plane. The relative deviation ^^7/7 of the width 7 
averaged over the trajectory is shown in grey scale (according 
to the colorbar). Additionally, the solid line indicates the bal- 
ance equation (|11|) and the broken line shows the equilibrium 
width of the ground state BEC in the given potential. 



3.2 Numerical results 

The above approach to the wave packet dynamics, which 
is based on the effective Hamiltonian Q, may be still 
oversimplified. For this reason and in order to check the 
estimate (|T^ we run the DNLS for different values of the 
nonlinearity g and harmonic confinement v. In order to 
reduce the number of independent parameters we have 
also assumed that the shape of the initial wave packet is 
defined by the ground state of the BEC before shifting 
the trap centre and, hence, the wave packet width 7 is 
no more an independent parameter (see the dashed line 
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in the stability diagrams in Fig. (S)). The results of these 
numerical studies can be summarised as follows: 

(i) The effect of stabilisation is sensitive to the initial shift 
of the packet Iq relative to the position of the separatrix 
0). In particular, no stabilisation was observed for Iq ^ 
/*. This is actually not surprising - the separatrix is 
the most fragile trajectory of the pendulum and any 
tiny perturbation destroys it. 

(ii) If Iq is sufficiently far away from the separatrix, there is 
a finite interval of nonlinearity ^min < 9 < ^max where 
the BEC oscillations are not decaying. 

(iii) The lower boundary ^min is defined by the condition 
for appearance of a (non-negligible) stability island for 
the effective system Q and is approximately given by 
Eq. (Pj) . 

(iv) The upper boundary ^max strongly depends on Iq and 
sometimes is not well defined in the sense that for a 
large g we find a transient or incomplete stabilisation. 
An example is seen in the lower right panel of Fig. ^ 
in the rapid increase of 1 — |?/^| for initial values near 
Iq ^80. This result (taken together with the existence 
of the stability island) suggests, that along with the 
stabilisation, the nonlinearity induces a different pro- 
cess in the system which destroys the regular oscilla- 
tions of the condensate when g exceeds some critical 
value. A more sophisticated approach than the vari- 
ational ansatz Q is required to take this effect into 
account. 

A boundary for the stability of dipole oscillations was 
proposed before in Ref. 5 by a simple argument based 
on the modulational instability of plane wave states. This 
argument lead to a critical value of l*/V2 for the ini- 
tial displacement. Our numerical calculations loosely sup- 
port this estimate for g/J^ 1 (see the upper right panel 
on Fig. EI) but also show additional dependence on g / J 
and r as well as significant deviations in other parameter 
regimes. It is important to realize, however, that instable 
motion below the separatrix (for /q < ^*) and the asso- 
ciated dephasing leads to a mean position of the wave 
packet at the equilibrium position of the harmonic po- 
tential whereas the stable or unstable motion above the 
separatrix (for Iq > /*) is characterised by a nonzero off- 
set from the equilibrium position. Hence both effects have 
very different character. 

3.3 Relation to Bloch oscillations of homogeneous 
BEC 

At this point we briefly mention the related problem of 
the dynamical (modulational) instability that has been 
studied in the context of the Bloch oscillations of a BEC 
subjected to a static force F. In this case the mean-field 
equation corresponding to Eq. {TJ reads 

ihai = (a^+i + ai-i) g\aifai + Fdlai , (13) 

where d is the lattice constant and the initial particle num- 
ber density |a^p = N/L is assumed to be constant and 




Ql , . . 1 

0.5 1 1.5 2 

J/Fd 

Fig. 6. Stability diagram for Bloch oscillations in a homo- 
geneous system. Here, ges = gN/L is an effective coupling 
constant. The dashed line separates the "universal" regime of 
weak forcing, where gcr — S.OFd, from the "non- universal" 
regime of strong forcing, where gcr additionally depends on J. 

large compared to one; N and L are the number of parti- 
cles and system size, respectively. Time-periodic solutions 
of Eq. correspond to Bloch oscillations and generalise 
the Bloch oscillations known from non-interacting parti- 
cles in periodic potentials. As shown in the recent papers 
|32ll33| . there is a critical value of the nonlinearity, be- 
low which these Bloch oscillations are stable, while above 
they are subject to a dynamical instability. This instabil- 
ity scrambles the relative phases and leads to inhomoge- 
neous, time-aperiodic density distributions. The stability 
diagram of Fig. El summarises the results. 

Since the dynamics of the atoms in a parabolic lattice 
for Iq ^ I* can be alternatively viewed as Bloch oscilla- 
tions in a static field with the local magnitude F = vl^/d^ 
a rough estimate for the expected instability regime may 
be drawn from this critical nonlinearity. However, an im- 
portant difference between the two systems is that the 
modulational instability analysis assumes a uniform state 
7 ^ oo, while in the pendulum dynamics the finiteness of 
7 is a crucial ingredient. Indeed, numerical explorations 
indicate that the parameter dependence is more compli- 
cated, which makes this an interesting problem for further 
study. 



3.4 Relation to superfluidity 

The concept of superfluidity rest on a rich phenomenology 
rather than precise definitions 34 . As mentioned above, 
we associate the emergence and breakdown of coherent 
dynamics with superfluidity in this paper in alignment 
with discussions in the recent literature related to cold 
atom experiments. In contrast to the traditional approach 
from condensed matter theory considering infinite systems 
in the thermodynamic limit, we are dealing here with an 
intrinsically finite dynamical system for which the concept 
of superfluidity yet has to be defined. 
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Clearly, a mere application of concepts borrowed from 
the theory of infinite systems will not help here: As an 
example we mention the Landau critical velocity, which 
is bounded from above by the speed of sound Vg, which 
is a function of the density. In a local-density-type argu- 
ment, we may consider the variation of the density over 
the wave packet and conclude that the critical velocity for 
the breakdown of superfluidity should vanish, as ^ 
in the tails of the wave packet. Thus, we would not ex- 
pect superfluid transport even though we implicitly as- 
sume that the atoms are completely Bose condensed by 
using Eq. 0. Nevertheless, we predict coherent transport 
in certain parameter regimes as discussed above. 

A systematic study of the robustness of the superfluid 
behaviour against energy dissipation from small impurities 
is beyond the scope of the current paper but will constitute 
an interesting extension of the present work. 

4 Conclusions 

In conclusion, we have considered the ID dynamics of a 
BEC of cold atoms in parabolic optical lattices. When 
interactions are absent, this system realizes the quantum 
pendulum Q with the experimentally controllable ef- 
fective Planck constant fieff = 2//*, where P of Eq. 
characterises the pendulum separatrix. The parameter P 
plays an important role both in theory and experiment. 
In particular, the relation between /* and the trap centre 
shift lo = Ax I used in the experiments to put the atoms 
in motion, defines whether BEC oscillations are symmet- 
ric with respect to the trap origin or not. The parameter 
fleff = 2/P also defines the rate of dephasing, because 
of which BEC oscillations decay even in the absence of 
atom-atom interactions. The effect of the latter on the 
discussed dynamics appears to be nontrivial. Naively, one 
would expect that any nonzero interaction enhances the 
decay of BEC oscillations. However, this is not the case - 
a moderate nonlinearity is found to stabilise the oscilla- 
tions, which now can be described in terms of the classical 
pendulum. The emergence of superfluid behaviour is thus 
related to a quantum-classical transition. We believe that 
for /o < ^* the stable regime of wave packet dynamics has 
been actually realized in the experiment 6 , where peri- 
odic oscillations of a BEC with a frequency given by the 
frequency of the classical pendulum have been observed 
and interpreted as a superfluid phenomenon. In order to 
see the transition to dephasing-dominated dynamics, the 
experiments would have to work at lower particle num- 
ber densities or reduce the nonlinear coupling constant, 
e.g. with '^Li atoms 35 , by tuning the atomic scattering 
length by means of a magnetic Feshbach resonance. Most 
surprisingly, stabilisation of wave packet motion may also 
occur above the separatrix, which appears to not have 
been observed yet. 
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